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Abstract 



C^ ' Exact reflection and transmission coefficients for supersymmetric sliape- 

^ • invariant potentials barriers are calculated by an analytical continuation of 

the asymptotic wave functions obtained via the introduction of new general- 
ized ladder operators. The general form of the wave function is obtained by 
/\ \ the use of the F(—oo, +00) -matrix formalism of Proman and Froman which 

C^ ' is related to the evolution of asymptotic wave function coefficients. 



Typeset using REVTgK 



*Electronic address: aleixo@nucth.physics.wisc.edu 
l^Permanent address . 

•^Electronic address: baha@nucth.physics.wisc.edu 
^Electronic address: macr@df.ibilce.unesp.br 

1 



I. INTRODUCTION 

Quantum tunneling through a potential barrier governs many interesting phenomena in 
physics ranging from fusion reactions in stars [Q] to the study of transitions from metastable 
states 0. There are very few exactly solvable examples of barrier penetration. Supersym- 
metric quantum mechanics has been shown to be a useful technique to explore exactly solv- 
able problems in quantum mechanics [Q. An integrability condition called shape-invariance 
was introduced by Gendenshtein ||^ and was cast into an algebraic form by Balantekin [Q. 
Reflection and transmission coefficients for a large class of shape-invariant potentials was 
given by Cooper et al 0. A general operator method for calculating scattering amplitudes 
for supersymmetric shape- invariant potentials was introduced by Khare and Sukhatme [|^. 
Even though an approximate method in the context of supersymmetric semiclassical approx- 
imation p] to calculate tunneling through one- dimensional potential barriers was presented 
in exact tunneling probabilities for shape-invariant barriers was not explicitly derived. 
We cover this latter subject in this article. 

Introducing the superpotential function 



W{x) 



h 



'im 






:i-i) 



where ^q[x) is the ground-state wave function of the Hamiltonian ff, and defining the 
operators 



A = W{x) + 



2m 



P, 



:i.2) 




it = w{x) - 

we can show that 

H-Eo = A^A . 
Since the ground-state wave function satisfies the condition 

i^o(a;) = 
the supersymmetric partner potentials 

Hi = A^A H2 = AA^ 



1.3) 



:i.4) 



:i.5) 



:i.6) 



have the same energy spectra except the ground state of Hi which has no corresponding 
state in the spectra of if 2- The corresponding potentials are given by 



Vi{x) = [W{x)f 



n dw 



;i-7) 



V,{x) = [W{x)f+ ^ ^^ 



(1.8) 



The shape-invariance condition [Q 

V2ix,ai) = Vi{x,a2) + Riai) (1.9) 

can also be written as |^ 

i(ai)i^(ai) = A\a2)A{a2) + R{ai) , (1.10) 

where ai 2 are a set of parameters that specify space-independent properties of the potentials 
(such as strength, range, and diffuseness) . The parameter a2 is a function of oi and the 
remainder R{ai) is independent of x and p. Not all exactly solvable potentials are shape- 
invariant [lin]. In the cases studied so far the parameters oi and 02 are either related by a 



translation [T3JTT[ or a scaling []12|. Introducing the similarity transformation that replaces 



Oi with a2 in a given operator 

f(ai)d(ai)ft(ai)=d(a2) (1.11) 



and the operators 



B+ = A\ai)f{ai) (1.12) 



B^ = Bl=f\ai)A{ai), (1.13) 

the Hamiltonian takes the form 

H-Eo = B+B^. (1.14) 

Using Eq. ( |1.10| ) one can easily prove the commutation relation 

[5_, fi+] = tt(ai)i?(ai)t(ai) = R{ao) , (1.15) 

where we used the identity 

R{an) = f{ai) Ria^^i) f\a,) , (1.16) 



valid for any n. Equation ( |1.15 ) suggests that 5+ and B^ are the appropriate creation 



and annihilation operators provided that their non-commutativity with R{ai) is taken into 
account. In this paper we extend the use of 5+ and B_ operators for calculating the 
asymptotic behaviors of the wave functions related with incidence of a particle on a super- 
symmetric shape-invariant potential barrier and obtain the exact transmission and reflection 
coefficients. 



II. EXACT WAVE FUNCTIONS 

The wave functions for all currently known supersymmetric shape-invariant potential 
barriers can be calculated analytically using supersymmetric operator techniques |TB|,|T3[. 
The final result can be expressed by single operators 5+ and I3z^ [0 or by couples of these 
operators. In the first case we can use the two additional commutation relations 



[B+B_,Bl] = Y.R{ak)Bl, (2.1) 



k=l 

and 



[fi+fi_,fin = 5;]i?(afc)Sl", (2.2) 



k=l 



obtained by induction using the relations 

R{ar,)B+ = B+R{ar,^^) , (2.3) 

R{an)B^ = B+R{an+i) , (2.4) 

that readily follow from Eqs. (|1 . 1 1|) , (|1.12|) and (|1.13|) . Considering that the Schrodinger 
equation can be write as 

fi+fi_^(x) =A^(x), (2.5) 

then the Eqs. ( p.l|) and (|2.2| ) imply that 5+ and 5_ can be used as ladder operators to 
solve the Eq. ( p. 51) P,0. To this end we introduce W_ \x) as the solution of the equation 

A_ (ai) ^L°^ (x) = = B^ (ai) ^L°^ (x) , (2.6) 

which implies 

,(0). . ( V2^ r 



vl/^^^(x,ai)~expl-^y d^Wi^,ai)\. (2.7) 

If the function 

n 

fin) = J2Ria,) (2.8) 

can be analytically continued so that the condition 

/(/x) = A (2.9) 

is satisfied for a particular (in general complex) value of fi, then the Eq. (|2.1|) implies that 
one possible form for the solution of Eq. (|2.5| ) is B^'^_ {x, ai) . Similarly if \E'^_ (x) satisfies 
the equation 

5+(ai)^f(x) = 0, (2.10) 

which implies that 

fia,)¥^\x) ~ exp h^fd^Wi^,a,)\ , (2.11) 

or 






;(0), 



(2.12) 



then the Eq. ( ^.21 ) imphes that other possible form for the solution of Eq. ( ^.51 ) is 
BZ^^ ^\_ (a;, ao) . At this point we conclude that the components of the wave functions, 
written in terms of the singles operators 5+ and Bz^ ■, for supersymmetric shape-invariant 
potential barriers can be written down as 



^+{x)=^BZ^-'^f{x,ao) 
where /i is obtained by the relation 



fc=i 



(2.13a) 
(2.13b) 



(2.14) 



where (3 and 7 are constants and 



vl/f (x,a^) =exp (±y^J^d^W{^,a,)\ . 



(2.15) 



With each value of fi we obtain several possible expressions for the components '^±{x) and 
the general expression for the wave function can be obtained with these components or a 
combination of them. 

It is also possible to express the components of the wave functions using couples of the 
operators 5+ and -Bz^. In this case we can use the equations (|2.1| ), ( p.2|) and the relations 
(|2.4D to show by induction that 



2n 



and 



[B+B^,iB^B 



\B+B^,{BZ'B, 



— iNn.l 



Y.R{ak){B+Bz'r^ 

k=l 



2n 



Y: R{ak) {BZ^B^r ■ 

k=l 



(2.16) 



(2.17) 



Using these last two equations and the same conditions ( |2.6| ) and ( p.lO| ) we can show 
that the components of the wave functions, written in terms of couples of the operators i?+ 
and Bz^, for supersymmetric shape-invariant potential barriers can be written down as 



^_(x) =(3{B+Bzy^'^^\x,ai) 
^+{x) =^{Bz'B^rBz'^f{x,ao) 

and where z/ is obtained by the relation 



2u 



A=Y^ R{ak) . 



(2.18a) 
(2.18b) 



(2.19) 



fc=i 



Note that in a given problem either Eqs. ( |2.13|) or Eqs. ( p.l8|) could be used, but not both 



III. ASYMPTOTIC WAVE FUNCTIONS 

The formal expressions for the components of the wave functions can be express in 
exphcit forms if we evaluate them asymptotically. In the case of single operators expressions 
we first note that using Eqs. (|1.12|) and ( |1.13|) the results (p.l3|) can be written as 



^_(x) = pA+{ai)A+{a2) ■ ■ ■ A+{a^)¥^\x,a^+i) (3.1a) 

^+(x) = -fAz\a,)Az\a2) ■ ■ ■ Az\a^+i)^f\x, a^+i) (3.1b) 

In this point we need to consider the two basic asymptotic behavior for the superpotential: 
i) W{x -^ ±00, a^) is constant (i.e., the potential barrier goes to a constant); and ii) 
W{x -^ ±00, a^) — > ±00 (i.e., the potential barrier goes to —00). In the former limit the 
commutator 



dx 



■^--,W{x,af,) 



W'{x,a^) (3.2) 

vanishes. In the ladder case this commutator can be ignored as 



/ W'(x a„ 

W{x,an)W{x,ak) + W'{x,an) = W{x,an)W{x,ak) [1+ ' 



W{x,an)W{x,ak)J 
^W{x,an)W{x,ak), (3.3) 

provided that W'{x, an) /W {x, an) remains finite, which is the case for all realistic superpo- 
tentials. Hence in both limits we can write Eqs. ( |3.1D as 

^_(x) = l3{Wi + W^+i){W2 + W^+i) ■■■{W^ + W,,+i) ¥'\x, a^+i) (3.4a) 

^+{x) = 7(1^1 + W^+,)-\W2 + W^+,)-' ■ ■ ■ {W,+, + W^+,)-' ^f (x, a^+i) . (3.4b) 

In these equations the quantity Wm is the short-hand notation for W{x, am)- 
If we assume that the superpotential satisfies the condition 

W{x,an) = W{x,ai) + {n-l)C{x), (3.5) 

then these asymptotic equations can be in a form suitable for analytic continuation, that is 

^-(^) =/5C'^ Y^^*^°^(^,«M+i) (3.6a) 

^^(x) = 7C-^-' T{2z + fi) .0)^ . ^3_g^^ 

+ ^ ^ '^ r(2z + 2/i + l) + ^ ' '^^' ^ ^ 

where z = W{x,ai)/({x). The condition given by Eq. (|3.5| ) is satisfied for a number of 
superpotentials and in the last section we give some examples. If this condition is not 
satisfied, the analytic continuation may still be done, but will be more complicated. 

When the asymptotic behavior of the superpotential is W{x, an) — > ±cx) we can use the 
identity 



to express Eqs. ( p.6|) in the simple form 



(3.7) 



(3.8a) 
(3.8b) 



We can repeat the same procedure used above in the case of couple of operators. Again, 
using Eqs. (|1.12|) and ( p..l3|) the results ( p.l8|) can be written as 



k=l 
u 

^+(a;) =7 n Az\a2k-i)A4a2k)Az\a2u+i)'^f\x,a2u+i). 



(3.9a) 
(3.9b) 



fc=i 



Considering the superpotential asymptotic simplifications given by Eqs. (|3.2| ), (|3.3| ) and 
the analytic continuation condition ( p.5| ) we can write the result for the components of the 
asymptotic wave functions in this case as 



^_(x) =p 



Til-z-u) 



,(0) 



T(l- z-2u)T(u + 



TT*^'(a;,a2.+i) 



,T, r^N . r(-z-2z/)r(z/ + |) (0) 



(3.10a) 
(3.10b) 



IV. GENERAL ASYMPTOTIC WAVE FUNCTIONS AND THE TRANSMISSION 

AND REFLECTION COEFFICIENTS 

Using the formalism developed in the Ref. |T^ we can write two possible asymptotic 
solutions for the one- dimensional time-independent Schrodinger equation in the form 



where 



with 



^i(x ^ ±oo) = An(±oo) /i(x) + A2i(±oo) /2(x) 

^2(X ^ ±00) = A2(±00) /i(x) + A22(±00) ^(x) 



/l(^) 



exp(-Nx(x)) 
q{x) 



and /2(3;) 



exp( -ix(x)) 
q{x) 



/x 
q{0 d^ and q{x) 

If we define the vectors 



'2m 
IT 



W{x,a^+i] 



(4.1a) 
(4.1b) 



(4.2) 



(4.3) 



*(x)) 






and 



then we can write 



(fix) 



(/iWI {f2ix)\] 



{^{x -^ ±00) 1= (fix) I A(±cx)) 
where the asymptotic coefficients matrix is given by 



v4ii(±00) Ai2(±oo) 

v42l(±00) A22(±00) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



Considering that the space evolution of the coefficients matrix A, given by an iteration 
process, can be written as 



A(x) = F(x,a;o)A(2;o) 



(4.8) 



then if we know the asymptotic coefficients in —00 and +cxd we can obtain the evolution 
matrix 



Using the Eq. 



F(-oo,+oo) = A(-oo) A ^(+00). 

we can show two basics properties of the F matrix: 

det F(— 00, +00) = 1 

F(-cx), +00) = F"^(+oo, -00) . 



(4.9) 



(4.10a) 
(4.10b) 



The knowledge of F(— 00, +00) permits the determination of the exact transmission and 
reflection barrier coefficients. If we consider a wave incidence from —00 to +00 then we can 
write the asymptotic wave function in the form 



or 









-^) = h{x) + Cnh{x) 
-00) = 6*7/2(2;), 



-00 J 
-00) 



(/(x) I a(-oo)) 
(/(x) I a(+oo)) 



where 



Considering that 



a(— 00)) 



1 
Cr 



and 



a(+oo)) 







(4.11a) 
(4.11b) 



(4.12a) 
(4.12b) 



(4.13) 



a(— 00)) = F(— 00, +00) \ a{+oo)) 



(4.14) 



we can conclude that 

Ct = TT-. r (4.15a) 

f..(-oo.+oo) 

Fi2(-oo,+oo) 

and the transmission and reflection coefficients are given by T =| Ct P and R =\ Cn p. At 
this point if we consider the conservation of probability, the time reversal invariance, and the 
invariance under space reflection it is possible to show the additional F-matrix properties: 

Fii(-oo, +00) = -F22(-oo, +00) (4.16a) 

F2i(-oo, +00) = -Fi2(-oo, +00) (4.16b) 

I Fi2(-cx),+cx)) I = I Fi2(+cx),-oo) |> 1 (4.16c) 

I F22(-oo,+cx)) I = I F22(+cx),-cx)) |<| Fi2(-oo,+oo) | . (4.16d) 



V. APPLICATIONS 
A. Parabolic Barrier 



For a parabolic potential barrier |[T8 



V^i(x)=K)-^mf]V, (5.1) 



the corresponding superpotential, obtained by using the Eq. ( |1.7|) , is given by 

W{x,ai) = aix (5.2) 

where Oi = ^iJ^Vt. The shape invariance condition ( [L.9[ ) imply that 

R{an) = ±2teo (5.3) 

where Eq = ^hQ. Using the Eq. ( p.l9|) we can conclude that 



J2 R{cik) = ±^2/i£o = A = E - K) T ^^0 (5.4) 

fe=i 



or 



1 A , , 

/x = --±«- (5.5) 

where A = (Vq — E)/eQ. The asymptotic form for the components of the wave function can 
be obtained using these results in the Eqs. 



^_(x)=/3(±i)-2±*2^±% 



exp (=F^g^) 
exp (±2f)^) 



(5.6a) 
(5.6b) 



or 



^„(x) = /3 



and 



^+(a;) = 7 



,T<fe+i)'r -(fc+i)7rA 



g±*(fc+f)^g(fc+|)^A 



,Ti(fc+i)7r (fc+i)7rA 



,±i(fc+f)^g-(fe+f)7rA 



f^^^t 



exp(=i=«f) 



^^*2 



v/^ 



exp(±2f)^ 



(5.7) 



\/^ 



(5.8) 



where q = JmVt/2hx and A; = 0, 1, 2, . . . . Using these results it is possible to write two 
asymptotic solutions of the Schrodinger equation when x -^ +oo in the form 



T / \ i^ -^ I iv^i exp(-«^^ 



v^ 



^2(3^ — * +CX)) = e *4e 4 I ^ 1*2 
therefore, if we identify 



exp(— «^ 



V^ 



+ e 4e 4 I ^ I '2 



+ e 4e" 



exp(+2^ ) 
exp(+i^^ 



/i(^) 



exp (—if) 



and /2(x) 



exp(+zf) 



we can conclude that the elements of A(+oo)-matrix will be 



74ii(+oo) = e 4e 4 | q 



pI 



yli2(+oo) = e *4e 4 I ^ 1*2 



In the case of x 



^2i(+oo) = e *4e 4 I ^ 
-00 if we consider that 



^22(+oo) =e*4e 4 I ^ I *2 



Q^'-^ = e 



=F(n+i)7rA I ^ |±2 
I c: I 



n = 0,l,2,..., 



(5.9a) 
(5.9b) 

(5.10) 

(5.11a) 
(5.11b) 
(5.11c) 
(5.11d) 

(5.12) 



in the Eqs. ( ^.8|) then we can write two asymptotic solutions of the Schrodinger equation 
when a; — > — cxD in the form 



^i(x 



-00 



+ 



e 4 + e 4)e 4 + e 4e 4 



■ TT . TT A ' TT TT A 

+ e 4)e 4 — e 4e 4 



^ 2 



^1 ^ 



exp(— zf) j 
exp(+i^^) 



+ 



(5.13) 



10 



and 



^lix 



-oo 



e 4 + e 4)e ^ + e ^e * 



g I 2 



exp[—ig 



+ 



[e 4 +e 4)e 4 



_,-7r SttA 

e 4e 4 



exp(+i^^) 



+ 



therefore we can identify the elements of A(— 00) -matrix as 



Ai(-oo) 

Al2(-00) 
^2l(-00) 
^22(-00) 



(e 4 + e 4)e 4 + e 4e 4 



e 4 + e*4)e 4 + e'^e 



3tt\ 

4 



e 4 + e 4)e 4 



e 4e 4 



e 4 + e 4 )e 4 — e 4e 4 



^ 



-t 



(5.14) 



(5.15a) 
(5.15b) 
(5.15c) 
(5.15d) 



In the choice of the two asymptotic wave functions for x -^ ±00 we considered the set of 
properties given by the Eqs. (|4.1CI|) and ([4.16|) that the F-matrix need to satisfies. Using the 
results for the A(— 00) and A(+oo) in the Eqs. ( |4.9| ) we can show that the evolution matrix 
can be written as 



F(— cx), +00) 



te 2 
-1 + te 2 ) 



Q 



-i\ 



'1 + ie 2 ) I ^ 

• 2LA 

le 2 



i\ 



and the exact transmission and reflection coefficients are given by 

1 1 



T 



and 



R 



'121 



'221 



-CX), +CX)) 



-CXD, +CX)) 



1 + e 



ttA 



,7rA 



'121 



-00, +00) 



l + e 



vrA 



(5.16) 



(5.17) 



(5.18) 



B. Morse Barrier 



For a Morse potential barrier |T9 

the corresponding superpotential, obtained by the Eq. (pTT]) , is given by 

W{x,ai) = ai + ae^/^ 
where 



ai = v^(l Tis) 



(5.19) 



(5.20) 



(5.21) 



11 



with e = h /{Smb"^) and s = JVq/s. The shape invariance condition ( |1.9| ) imply that 



where a„_|_i = a„ + 2^/e. Using the Eq. (p.l9| ) we can conclude that 

Y^ R{ak) =a\- a^+i = A = £■ + a^ 



fc=i 



or 



If we remember that 



a.+i = ±.i\fE . 



/i 



Cf(t+1 ~ '^1 



2v^ 



we can use the Eqs. ( ^.21| ) and (|5.24]) to show that 



1 s r 
u = — ±%- ±%- 
^222 



(5.22) 



(5.23) 



(5.24) 



(5.25) 



(5.26) 



where r = J E /e. 

Considering the asymmetry of the Morse potential barrier the wave function will have 
a different behavior in +oo and — oo. Therefore the asymptotic form of the components of 
the wave function for x — > +cx3 can be obtained using the last results in the Eqs. (p.8[) 






exp (^=Fif exp(a;/6)j 



exp(s/6) 
exp (±«f exp(x/6) 
exp(x/6) 



(5.27a) 
(5.27b) 



Using these results it is possible to write two asymptotic solutions of the Schrodinger equa- 
tion when X -^ +oo in the form 



^i(x -> +oo) = e^^e'^^''^''^ 



exp (—if exp(x/6) 



exp(2;/6) 



(5.28a) 



+ e 4e 



i-..-As-r) 



exp (+i|exp(a;/6) 



exp(a;/6) 



^2(3: -^ +00) = e"*5e"?(''"''^ 



exp (—if exp(x/6) j 
exp(x/6) 

exp f+if exp(x/6) 



(5.28b) 



+ 6*6 



ifpf(«+0 



exp(x/6) 



12 



therefore, if we identify 

exp(-i|exp(x/6) 
/i(x-^+oo) = — 



exp(x/6) 



and /2(x -^ +00) 



exp (+«f exp(x/6) 



exp(x/6) 



e-fe-f(-0 



we can conclude that the elements of A(+cx3)-matrix will be 

Au{+oo) = e^f e-5(^+") 

^i2(+oo) 

A2i(+oo) = e-'?e?(^-") 

A22{+oo) = e'^e^^'+'K 

In the case of x ^ —00 we can substitute for s and r and using the Eq. (p.6|) find 



^-(x 



-oo) = /3^^^I^i^^e^- 



\l'+(x — > —00) = 7 



r(i±z|± 



i/r r(±2r 



2/ ±ifcx 



(5.29) 



(5.30a) 
(5.30b) 
(5.30c) 
(5.30d) 



(5.31a) 



(5.31b) 



where k = y2rnE/h. Using these results it is possible to write two asymptotic solutions of 
the Schrodinger equation when x — > — 00 in the form 



^i(x 



-00 




(5.32) 



and 



^2(x 



-00 



e 4 



+ 



6 4 62 



(s+r) 



rQ + ^f + ^i) r(i + ^f 



+ 



therefore, if we identify 



6 4 6 2 



+ 



e' 4 e 2 



r(l-^f-zi) r 



^2 ^2 



/i(x ^ -00) 



exp (+iA;x) 
7k 



and /2(x ^ —00) 



exp {—ikx) 



(5.33) 



v^ 



(5.34) 



we can conclude that the elements of A(— oo)-matrix will be 



13 



An(-oo) 
^i2(-oo) 

^2l(-00) 
^22(-00) 





e 


' 4 e 2 




r 


(^ 


+ if + 


'§) 


r 


(^ 

e* 




) 


r 


(^ 


■ TT Trr 

g*4 e 2 


■i) 



ra 



^2 



+ 



+ 



+ 



+ 



6 4 e 2 



6 4 62 



— «7 



f(^+r) 



6 4 



V2 ^2 "'" ^2 

■ Stt 7rs_ 
6 4 6 2 



ra 



^2 



■\/r r (— ir) 
i/r r (ir) 
■\/r r (— ir) 



Again, in the choice of the two asymptotic wave functions for x -^ ±oo we have 
ered the properties of the F- matrix. Using the results for the A(— oo) and A(+cxd 
Eqs. ( |4.9| ) we can show that the evolution matrix can be written as 



F(— oo,+oo) 



ig ih 
ih* ig* 



where 



6f(^-^)^r(i 



tr] 



and 



r(l + ^f + ^i 



6f(^+'-)^r(2r) 



r(i + 



''2 2 



In this case the exact transmission and reflection coefficients are given by 

1 6"^'-''"'"''^ sinh (vrr) 



and 



R 



1 Fi2(-00,+CX)) |2 


cosh 


'_^{s-r) 


; 


F22(-00,+00) 2 


6"^^^ cosh 


f(s+^y 


i^l2(-00,+00) |2 


cosh 


TV 

2 


(s-r) 





(5.35a) 
(5.35b) 
(5.35c) 

(5.35d) 

consid- 
in the 

(5.36) 
(5.37) 



(5.38) 



(5.39) 



(5.40) 



C. Eckart Barrier 



For an Eckart potential barrier |^ 



Vi{x) = Vosech^{x/2b), 
the correspondent superpotential, obtained by the Eq. ([U^ ), is given by 



(5.41) 
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W{x, ai) = ai tanh {x/2b) (5.42) 

where 

ai = v^(-l±is) , (5.43) 



with e = h /{32mb ) and s = JVo/e — 1. The shape invariance condition ( pToD imply that 

R{an) =al- a^+i , (5.44) 

where a„+i = a„ — 2^/e. Using the Eq. (|2.19|) we can conclude that 



2u 

J2 R{ak) =al- a^+i = A = E + aj (5.45) 

fc=i 

or 

a2.+i = ±2\/^. (5.46) 

If we remember that 

2.= ^^^i^, (5.47) 



we can use the Eqs. ( ^.43| ) and (|5.46|) to show that 



2u = -^±t'-Ti^^, (5.48) 

where r = J E/e. 

The asymptotic form of the components of the wave function can be obtained using these 
results in the Eqs. ( p.lO| ) 

^_(x) = 13 —, ^ J ^ — ^ e^*^" (5.49a) 



r( 


;!±^t±^i) 




r(i±i 
r(±^i) 







vl/+(x) = 7 ^ JL , \: , \ ^' e^''- , (5.49b) 



r(!±^i± 



where k = y/2mE/h. Using these results and the relation 

v4 / v4 / cosh (vry) ± sinh (vry) ' 

it is possible to write two asymptotic solutions of the Schrodinger equation when x — > +cxo 
in the form 

^i(a; ^ +oo) = Cie-^'^^ + C*e^'=^ (5.51a) 

^2(3; -^ +00) = C*e-*^^' - Cae*'^^ , (5.51b) 
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where 



Ci = v2 TT < cosh 



- (r ± s) 



isinh 



-(r±s) 



-1 



xir|i + .- r -±,--.j r -T.j-.j 



-1 



and 



C2 = V2nT [i-] <!cosh 



fr ± si 



, ^ /3 s r\ ^ /3 s 
X Jr-^z- + z- r-±z- + « 
^ 14^4 4/ V4 4 4 



isinh 
r 



TC 



(r ± si 



Therefore, if we identify in the Eq. ( 5.51 

„—ikx 



\-ikx 



fi{x -^ +oo) 



Vk 



and f2{x -^ +oo) 



y/k 



we can conclude that the elements of A(+oo)-matrix will be 

An{+oo) = Ci 

Aui+oo) = c; 

A2i(+oo) = CI 



(5.52) 



(5.53) 



(5.54) 



(5.55a) 
(5.55b) 
(5.55c) 
(5.55d) 



Also we can write two asymptotic solutions of the Schrodinger when x — i> — oo in the form 



^i(x 
*2(a; 



-oo) = Cie'^'' - C*e 

-OO; = — OoC — 02 



* „—ikx 



C* ikx /^ —ikx 
2^ — <-^2e , 



and identifying 



fiix 



-oo 



+ikx 



Vk 



and /2(x — > — oo) 



-ikx 



Vfc 



(5.56a) 
(5.56b) 



(5.57) 



we can conclude that the elements of A(— oo)-matrix will be 



yln(-oo) = Ci 

Ai2(-oo) = -c; 

A2i(-00) = -CI 
^22(-00) = -C2 



(5.58a) 
(5.58b) 
(5.58c) 
(5.58d) 



Using the results for the A(— 00) and A(+oo) in the Eq. ( [4. 91) we can show that the evolution 
matrix can be written as 



F(— 00, +00) 



9 

-h* 



h 
-9* 



(5.59) 
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where 



and 



h 



O1U2 — L'iL'2 
C1C2 + CxC2 

2C1C2* 



(5.60) 



(5.61) 



0102 + '-^l*-^2 

On using the Eqs. (|5.52|) , (|5.53|) after a considerable amount of algebra we can show that 
the exact transmission and reflection coefficients are given by 

1 _l Ch Cl^ sinh^ (f ) 

4 ci^'cl 



T 



Fi2(-oo,+oo) P 



and 



R 



'221 



-CXD, +00) 



'121 



-00, 



-CX3 



9i 






sinh^ (^f) + cosh^ (^ 

cosh^ (f ) 
sinh^ (f ) + cosh^ (f ) ■ 



(5.62) 



(5.63) 



VI. CONCLUDING REMARKS 

In conclusion we note that the present technique is a powerful and an elegant prescription 
to obtain exact reflection and transmission coefficients. This method may also be used for 
all supersymmetric shape invariant potential barriers that satisfy the analytic continuation 
condition ( |3.5| ). 

One possible application of shape-invariance formalism is to multidimensional quantum 
tunneling. In nuclear physics applications multidimensional quantum tunneling can be vi- 
sualized as the tunneling of a quantum mechanical system (such as a nucleus with internal 
excitation) instead of a structureless particle through a one-dimensional barrier. The nucleus 
typically taken to enter the barrier in its ground state and may emerge either in the ground 
state or in an excited state at the other side of the barrier. The interaction between the 
penetrating quantum system and the barrier also needs to be specified based on the physical 
conditions of the problem. It has been known for some time that, if the excitation energies 
are neglected, the penetration probability of an A^- dimensional system can be reduced to a 
sum of probabilities of A^ one-dimensional suitably defined barriers [^. The eigenchannel 
formulation remains valid even for finite excitation energies as long as the energy-dependence 
of the weight factors is taken into account p2[. Our formulation would be applicable in such 



cases if the eigenpotentials are shape- invariant. A simpler limit would assume factorization 
of the interaction between the barrier and the quantum system into a product of two quan- 
tities which are functions of the barrier and internal degrees of freedom respectively. Such 
a factorization approach was already applied to a coupled system of equations for bound 
states 



Our formulation also casts the tunneling problem in an algebraic basis 0,0. If the 

internal system can be described by an algebraic model such as the interacting boson model 

23| then it may be possible to cast the entire problem into an algebraic framework. A 

group-theoretical formulation can be a starting point of systematic approximations such as 



those given in Ref . [^ . A detailed study of such aspects is deferred to later work. 
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